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Abstract. For a general class of linear collisional kinetic models in the torus, 
including in particular the linearized Boltzmann equation for hard spheres, the 
linearized Landau equation with hard and moderately soft potentials and the 
semi-classical linearized fermionic and bosonic relaxation models, we prove ex- 
plicit coercivity estimates on the associated intcgro-diffcrential operator for some 
modified Sobolev norms. We deduce existence of classical solutions near equilib- 
rium for the full non-linear models associated, with explicit regularity bounds, 
and we obtain explicit estimates on the rate of exponential convergence towards 
equilibrium in this perturbative setting. The proof are based on a linear energy 
method which combines the coercivity property of the collision operator in the 
velocity space with transport effects, in order to deduce coercivity estimates in 
the whole phase space. 
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1. Introduction 

In this paper we study a general class of linear inhomogeneous kinetic equations 
in the torus (including linearized Boltzmann, Landau, classical relaxation, semi- 
classical relaxation and Fokker-Planck equations). We then use the properties ob- 
tained on their evolution semi-groups to gain insight into the behaviour of the full 
non-linear models in the cases of Boltzmann, Landau and semi-classical relaxation 
equations. The main tool is an estimate of the following type. Assume that the 
linear (self-adjoint) collision operator L acting only on the velocity space is coercive 
for a certain norm, and that L has a structure "mixing part + coercive part" (which 
shall be given a precise meaning in Subsection I1.2J1 . Then the integro-differential 
operator T = L — v ■ V x . (taking into account transport effects) acting on the phase 
space of positions x and velocities v, satisfies a coercivity property in x and v, which 
implies in particular the existence of a spectral gap estimate when L has a spectral 
gap in velocity (note that in general T is not sectorial). Moreover our proof shows 
how to compute the constant of coercivity of T according to the one of L in the 
velocity space (and thus the rate of exponential convergence to equilibrium when L 
has a spectral gap in velocity). Before we explain our method and results in more 
detail, let us introduce the models and problems in a precise way. 

1.1. The Problem and its motivation. In order to study the convergence to equi- 
librium, a new quantitative method in the large, the "entropy-entropy-production" 
method (EEP-method), has developed from the beginning of the 1990 decade, see 
EH] in the spatially homogeneous setting, and [TT] in the spatially inhomo- 
geneous setting. It has provided new powerful and robust tools in the study of 
relaxation towards equilibrium (see for instance jSlHH El [13 EH] ), and seems to be 
the best suitable approach to deal with non-linear models in the large. However, it 
can be seen from these references that the method, while very robust with respect 
to nonlinearities and able to deal with external potentials at no increased difficul- 
ties, has two major shortcomings. First it relies on uniform in time estimates of 
regularity of the solution, that are usually hard to establish. Second, it seems to fail 
to give the optimal rate of convergence, in particular when L has a spectral gap in 
velocity it seems to fail to give exponential rate of convergence to equilibrium. 

In this article we provide a linear energy method in order to overcome these prob- 
lems for linear models, or for non-linear models in the perturbative setting (this 
method is explicit in terms of the explicit coercivity estimates on the linearized col- 
lision operators in the velocity space, see P3 E21 EII| ) • ^ ne difference when compared 
to the EEP-method is that the convergence as well as the uniform in time regularity 
bounds are obtained in a single step. However this approach is linear and therefore 
limited to the perturbative setting near equilibrium for the non-linear models. A 
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similar approach has previously been used in Guo's papers [2U1 El 1221 EZ|- Also 
another approach has been developed in [2H] for the Fokker-Planck model in a con- 
fining external field (see also [21]), and Herau recently generalized the method to 
the linear relaxation model in j2H]- Similar results for the Fokker-Planck equation 
were recovered in jH] by an energy method. The results we present here are partly 
included in these references, however the proof is simpler and more explicit. Our 
viewpoint unifies previous scattered results and explores new situations such as semi- 
classical relaxation models. The works quoted so far are the starting point of this 
paper. In particular one of the key ideas of our proof, namely looking for the time 
derivative of a "mixed term" of the form J V x h ■ V v h, is inspired from 03] where it 
first appeared in this explicit form (see also 

Most of these articles deal with the Boltzmann equation that we address in sec- 
tion 15.41 Our main abstract theorem applies to the Boltzmann linearized colli- 
sion operator for hard spheres (and hard potentials with Grad's cutoff assumption): 
smooth solutions with explicit regularity bounds and rate of convergence to equi- 
librium are constructed near equilibrium. A similar study has been performed in a 
non-constructive way for the Vlasov-Poisson-Boltzmann system near equilibrium in 
[TS] and for the Landau equation in ^U]- I n [2H] the same method was also applied 
to the Boltzmann equation for cutoff soft potentials. 

Guo's argument relies, roughly speaking, on the coercivity of the linearized Boltz- 
mann operator "in the mean" - i.e., if integrated over a time interval - for small 
perturbations of the Maxwellian. Due to this averaging in time precise rates of con- 
vergence are quite hard to deduce from the results in ^SJ|TI|1 . Later Guo refined his 
approach in [21] to meet the needs of application to the Vlasov-Maxwell-Boltzmann 
system. The new idea is to use a norm that includes temporal derivatives. In this 
norm the fields are almost controlled by the deviation of the distribution from the 
Maxwellian - to be more precise they "loose" one derivative. This leads to instan- 
taneous coercivity and thus global perturbative solutions for the Vlasov-Maxwell- 
Boltzmann system as well as exponential decay to equilibrium in the case of Vlasov- 
Poisson-Boltzmann, in these norms. The proof, while being very complicated, is 
constructive. More recently a farther refined method based also on norms including 
temporal derivatives has been used to show almost exponential decay to equilibrium 
for various of the mentioned models and the relativistic Landau Maxwell system in 

The EEP-method has been applied to the Boltzmann equation in ^2] ~~ again 
assuming uniform in time regularity bounds on the solution - but for a large class 
of collision operators and most important - without perturbative assumptions. All 
these works have been carried out on bounded domains, essentially the torus apart 
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from ^2] where various types of boundary conditions are considered. For complete- 
ness we also refer to a recent preprint dealing with the problem in whole space j22j 
and the article [2Zj for a related non-linear energy method. 

Despite this vast amount of recent literature on the problem our proof has two 
advantages. First it is simpler than the ones in the articles quoted above, mostly 
due the fact that we take advantage of the mixing properties of the collision oper- 
ator. Second it separates in a very clear way between linear effects in these energy 
estimates (transport + linear collision), which are expressed in a coercivity estimate 
on this linear part, and the problems arising from the small remaining bilinear part 
when considering solutions near the equilibrium. We are able to derive exponential 
convergence to equilibrium without resorting to norms including time derivatives 
and in a purely instantaneous manner. 

1.2. The models. We will study initial value problems for equations of the form 

(1.1) dtf + v vj = Q(f), t > o, x g n, v G R N , 

where N > 1 denotes the dimension of space. In this equation, / denotes the 
distribution of particles in phase space, therefore it is a time-dependent non-negative 
L 1 (fi x ~R N ) function. The operator Q models the collisional interactions of particles 
(either binary or between particles and a surrounding medium). It is local in x 
and t and it depends on the particular model of interaction chosen. The term 
v ■ V x / corresponds to the free flow of particles. For the spatial domain Q, we shall 
consider here the periodic case, that is Q = T N . Hence / = f(t, x, v) satisfies initial 
conditions f(t = 0,x,v) = f (x,v) and periodic boundary conditions. 
One defines the global mass, momentum and energy of the solution / as 

/ fdxdv, / fvdxdv, / / |t>| 2 da; dt>. 
Jqxr n JnxR N JnxK N 

Depending on the model chosen for the collision operator Q, one has preservation 
along time of part or all of these quantities. 

Under very general assumptions, equations of class (|1-H) admit a unique equilib- 
rium in the torus, which we shall denote by f^, and which is independent of t,x 
(this is trivial for relaxation models admitting only mass conservation, and for mod- 
els admitting conservation of mass, momentum and energy, such as Boltzmann and 
Landau equations, this is shown easily inspiring from the arguments in ^2] and 

Then one can consider the linearization around this equilibrium. In order to 
reduce to a Hilbert space setting, one usually considers perturbations of the form 

1/2 

/ = foe + foa h. Discarding the bilinear term, it yields the following linearized 
equation on h 

(1.2) d t h + v-V x h = L(h), 
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where L depends on the precise form of the collision operator Q. The unknown h 
belongs to L 2 (T N x R N ). 

Now let us give a general framework for linear collisional kinetic models. We shall 
denote as usual by L 2 the Lebesgue space of square integrable functions, and, for 
k G N, H k the Sobolev space of L 2 functions with square integrable derivatives up 
to order k. When needed we shall indicate as a subscript the variables (x or v) these 
functional spaces refer to. For the variable x, Sobolev spaces refer implicitly to 
periodic Sobolev spaces. When no subscript is used, these functional spaces always 
refer to x and v. 

HI. ( Structure ) We consider a linear collision operator L on L 2 = L 2 (T N x M. N ) 
which a closed and self-adjoint operator on L 2 and local in t, x. 
We assume that it writes 

(1.3) L = K — A 

where A is a coercive operator in the sense: there is a norm || ■ ||a„ on M. N 
(the space of velocities), such that 

\\h\\ 2 L 2 < Vi\\h\\ 2 Av < (A(h), h) L 2 < ||/i||a„ 

and also 

(V v A(h),V v h) Li > ^3 A ||V>||l-^ A |l^lli ? 

for some constants v§ , v^, u^, > 0. 

Moreover we assume that L satisfies (for some constant C L > 0) 

(L{h),g) Ll <C L \\h\\ K \\g\\ K . 

To shorten the notation we also introduce the norm 

H2. (Mixing property in velocity) We assume that K has a regularizing effect 
in the following sense: for any 5 > 0, there is some explicit C(S) such that 
for any h G H^, 

(1.4) (V v K(h),V v h) L2v <C(5)\\h\\l i + 5\\V v h\\h ■ 

H3. (Relaxation to the local equilibrium) We assume that L, as an operator on 
Ll, has a finite dimensional kernel 

N(L) = Span {</?i, . . .,</?„} 
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and we denote by Hi the orthogonal projection on N(L) in l? v . We make the 
following local coercivity assumption: there is A > such that 

(L(h), h) L 2 < -A \\h - Ui(h)\\\ v . 

This together with HI implies in particular that L is non-positive and has 
a spectral gap in L 2 , whose size is bounded from below by 

a 1= (*)a>o. 

In the structure assumption HI, K shall typically stand for a multiplicative oper- 
ator by a function v (usually called the collision frequency) in the case of short-range 
interactions, or some diffusion operator in the case of long-range interactions. The 
norm || • ||a that we define can also loosely be seen as the norm of the graph of A 1 / 2 . 
The coercivity of A is linked with the coercivity property of the whole linearized 
collision operator L in H3 (where the word "local" refers to the position x and the 
fact that it pushes the dynamic towards local equilibrium). This coercivity property 
is crucial to ensure exponential decay towards equilibrium, even in the homoge- 
neous setting. In the case where it is weakened, which happens for interactions with 
"weak collision effect" such as soft potentials (see for instance in |H2])> one expects 
convergence rates of the form e~* T with r < 1 (see the results in |38j ) . 

The kernel of L in L 2 , which is composed of functions belonging to N(L) for any 
x, corresponds to the manifold of local equilibria for the linearized kinetic models. 
Therefore when the x variable is added, II; is the projection on the "fluid part", and 
(Id — II;) is the projection on the "kinetic" part. It is defined by 

n 



i=i 

It is trivial that any local equilibrium uniform in space is indeed a global equilib- 
rium. Since L is self-adjoint, the (fx, . . . , <p n belong to the kernel of its adjoint L*, 
and thus by integrating in x and v we get 

d f 

Vi = l,...,n, — / htpi dxdv = 0. 



dt 



Hence we denote 



n g (M = ( / hipidxdv) 

~l \Jt n xR n J 



which is time and space independent (this can be shown easily to be the orthogonal 
projection on N{T) in L\ v , directly or using the coercivity property of Theorem ll.lj) . 
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A detailed study of physical models satisfying assumptions H1-H2-H3 shall be 
given in Section 

1.3. Main results. First we state the main result on the coervivity estimates on 
T and the consequence on its evolution semi-group. 

Theorem 1.1. Let L be a linear operator on L 2 XV satisfying assumptions H1-H2- 
H3. Then T = L — v ■ V x generates a strongly continuous evolution semi-group e tT 
on H\ v , which satisfies 

||e tT (Id-n g )||^ <C T exp [-vrt] 

for some explicit constants Ct, tt > depending only on the constants appearing 
in assumptions H1-H2-H3. More precisely 

VheH 1 , (Th, h) w < -C' T (\\h - U g (h)\\l + \\V x , v (h - U g (h))\\l) 

for some (explicit) Hilbert norm 7i l equivalent to the H 1 norm, and some explicit 
constant C' T > 0. 

Remarks: 

1. Note that under slightly strengthened assumptions, a similar result is proved 
in H k in Theorem 13. 11 

2. The method does not rely on an abstract result from spectral theory such as 
Weyl's theorem, like Ukai's proof of the existence of a spectral gap for the Boltzmann 
equation for hard spheres in [1UJ. Hence we do not need the compactness property 
of K, although we require a regularizing property on K which is strongly related 
(see the discussion in [HJ). Our method can be seen as a quantitative version of 
Ukai's result (in the case of the linearized Boltzmann equation). In particular it 
shows that apart from 0, the spectrum of T is included in 

{£eC;Re(0<-<r T }. 

The abstract setting emphasizes what is effectively required from the linearized 
collision operator to deduce exponential convergence, and it allows to apply the 
method to other models as well. Since for the linearized Boltzmann equation it was 
proved in [30] that T is not sectorial (its essential spectrum is given by a half-plane), 
our work can also be seen as a method to prove exponential decay of the semi-group 
for a class of non-sectorial operators, which is in general quite tricky. In the case 
of Fokker-Planck type operators, other methods have been developed in (23 Hlj to 
solve this question. 

3. In order to obtain a completely quantitative result of convergence to equilib- 
rium, one has to get estimates on the constant A in assumption H3. This ques- 
tion had remained open for a long time for important physical models such as 
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the linearized Boltzmann collision for hard spheres or the Landau collision oper- 
ators for hard and moderately soft potentials. It has been solved recently in the 
works m E2 IHHj- Therefore this theorem allows to compute rates of convergence 
explicitly for all the models we consider in Section |S] (except semi-classical relaxation 
for bosons). 

In the two following theorems we define ko the smallest integer such that E{ko/2) > 
N/2 (where E denotes the integer value). Our second main result is for the nonlinear 
Boltzmann and Landau models. 

Theorem 1.2. Let us consider either the Boltzmann equation ()5.7|) for hard spheres 
or hard potentials with cutoff' or the Landau equation (j5.19|) with 7 > —2, in the 
torus. Let < fo be an initial datum with finite mass and energy, and we denote 
by foo the unique equilibrium associated to fo. We suppose that the initial datum 
satisfy s 

ll/oo 1/2 (/0-/co)||^<£ 

for some k > k and some < e < £0 where eo depends explicitly on the collision 
operator. 

Then there exists a unique global solution < / = f(t,x,v) G C([0, oof, H k ) of 
the initial value problem (jl.ljl . such that 

Vi > 0, ||/- 1/2 (/(t, •, •) - /oo)|U < C exp l-rt] 

for some explicit constants C, r > 0. 

The conclusion still holds true when a repulsive self- consistent Poisson potential 
is added (without smallness condition on the intensity of the self- consistent interac- 
tion), in the case of the Boltzmann equation for hard spheres, or the Landau equation 
with 7 > — 1. 

The next theorem deals with the semi-classical relaxation models. 

Theorem 1.3. Consider the semi- classical relaxation equation (J5.2)) in the torus 
for fermions (e = 1) or bosons (e = —1), with an initial datum < fo G L 1 . Let the 
equilibrium distribution be given by 

1 + e /too M 

with 1 + eKooAi > 0, where M. is the normalized Maxwellian 

P -M 2 /2 

M(v) 



(2tt) n / 2 

and is defined by mass conservation. Let k > ko and let the initial datum satisfy 
(l + eK 00 M)(K oc My 1/2 (fo-fo 



< E 
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for some < e < Eq where Eq depends explicitly on the collision operator. 

Then there exists a unique global solution < / = f(t,x,v) G C([0, oof, H k ) of 
the initial value problem (jl.ljl . such that 



Vt > 0, 



1 + e Koo M) ( Koo A^) _1/2 (/(*, •, •)-/») < C exp [-rt] 



/or some explicit constants C,t > 0. 



Remark: The condition on the form of the equilibrium distribution is in fact 
trivially fulfilled for the fermionic case. In the bosonic case it is a condition of 
smallness on and thus on the initial mass. Indeed it is equivalent to impose that 
the mass of the initial datum is small enough such that no condensation occurs. 
This is not for technical reasons but necessary to ensure exponential convergence as 
can be seen from the detailed study of the asymptotics in the spatially homogeneous 
case in [T3] . 



1.4. Outline of the article. The article is structured as follows. In Section |2] we 
give the proof of Theorem ll.il which turns out to be very short and simple. Then in 
Section El we expose several extensions of the method. In particular we show how to 
generalize Theorem 11.11 to higher-order Sobolev norms, and how to include a weak 
external field or a self-consistent Poisson potential in the study. Section 0] is devoted 
to the application of the previous study to the genuine non-linear problems of the 
form (jl.ljl near equilibrium, and we prove the abstract Theorem 14.11 Finally in Sec- 
tion we discuss the general assumptions of Theorem 11.11 and Theorem 14. II for an 
extensive list of physical models: classical or semi-classical relaxation, Boltzmann 
equation for hard spheres or hard potentials with cutoff, Landau equation for hard 
potentials or moderately soft potentials, linear Fokker-Planck equation. Then The- 
orem 11.21 follows from this study together with the Theorem 14.11 (we also comment 
on the marginal differences between the proofs of Theorem 11.21 and Theorem I1.3JI . 



2. Proof of Theorem 11.11 

We divide the proof into several steps. 

Step 1. Since by assumption H3, the operator L is non-positive, the L xv norm 
is decreasing along the flow and it is straightforward to deduce that T generates a 
strongly continuous contraction evolution semi-group on L 2 XV (see for instance). 
In order to estimate the semi-group in H xv) let us consider h G H xv n Dom(T) fl 
Ni^T) 1 - and h = h(t,x,v) the associated solution of the equation dth = T{h). Now 
we study the evolution of the H]. v norm of h. 
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Step 2. We estimate the time evolution of the L 2 norm of h. Using the skew 
symmetry of the transport part together with assumption H3, we find immediately 



(2.1) 



d 
dt 



Step 3. We estimate the time derivative of the gradients in x and v. 

• For the gradient with respect to x we obtain, thanks to assumption H3, 



(2.2) 



^||v^ni 2 < -2A||v x /i-n,(v x /oilA- 



For the gradient with respect to v we get 
d o 

— || V v h\\l* = 2 (V v K(h), V v h) L 2 - 2 (V v A(h), V v h) L 2 - 2 (V x h, V v h) L 2. 



Using assumption HI we have 



-2(V v K{h),V v h) L 2 < -2z/ 3 A ||V^||i + 2z/i 



2 

L 2 ' 



Assumption H2 yields 



z/ A z/ A 



2 (V v K(h), V v h) L 2 < -L _L \\V v h\\l 2 + 2 C 



±± 
4 4 



2 < 

L 2 - 



|V>||i + 2C 



±4 



2 

L 2 



Furthermore we split 



2 (V x h,V v h) L 2 < V -\\ \\W v h\\l 2 + ^4 WVxhWl, 



A 2 



^ 2 



K ^3 



Using the last three inequalities we have 
Now we write 



<fllv,* + ^l|v,C 2 



i 2 + jSll v ^ll' 2 - z/ 3 A llv>|iL 

^1 ^3 



l/ll 



L 2 



< 2 ll/i-IIi 



2 2 + 2 wmww* ■ 



Since IT s (/i) = we deduce that Tli(h) has zero mean on the torus, and Poincare's 
inequality in the torus yields (for a constant Cp only depending on the dimension 
N) 

\\^i{h)f L 2 < CpWU^VMW 2 < C P \\V x hf», 
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and thus we get for some explicit constants Ci, C 2 > 

(2.3) ^llV^Ula < C x \\h - II,(/0|ft + C 2 \\V x h\\h - 4 \Wvh\\l 

• For the mixed term we have 

^(V x h,V v h)& = -\\V x h\\ 2 & +2(V x L(h),V v h) L 2. 
Then we write (using assumption HI and Cauchy-Schwarz's inequality in x) 

2(V x L(h),V v h) L 2 = 2(L (V^-rLXV^))^/*)^ 

< 2C L nv^-n^v^)^ \\v v h\\ A 

< C L V \\V x h-Tl l (V x h)\\l + C L T ] - 1 \\V v h\\l 

for any i] > 0. Hence we obtain 



^<VA V v h) < - \\V x h\\ 2 L2 +C L V \\V x h - Ik(y*h)\\l + C L rT x \\V v h\\ 2 A . 



Step 4. Now it remains to combine equations ()2.1|) . (|2.2|) . ()2.3|) and (|2.4|) : we pick 

A, a, /3, 7 > and compute 

(l [ A\\hf L2 + a\\V x h\\ 2 L 2+ (3\\V v h\\h+l(^ x h,V v h) L 2 



dt 



< (PC, - 2AX) \\h - H^Wl + ( V1 C L - 2a\) \\V x h - Ui(v x h)\\l 



+ {q- l iC L - ||v^||i + (c 2 /5- 7 ) nv^iila. 

For a given (3 with u^(3 > 1, we first fix A big enough such that 

{(3d -2 AX) < -1, 

then 7 big enough such that 

(c 2 /?- T )<-i, 

then 77 big enough such that 

(r l ic L - (3i4) < -i, 

then a big enough such that 

(r]jC L - 2aA) < -1 
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and such that 7 2 < a/3 and a > (3. For this choice we obtain 
d 



dt 



A \\h\\ 2 L2 + a \\V x h\\ L i + P \\V v h\\ 2 L i + 7 (V s /i, V v h) 



< 

The function 



v x h\\h + liv^lli + \\h - n,(/i)||i + ||v x /i - n,(v^)|| A 



.F(i) = [A ||/i||i 2 + a \\V x h\\h + P II V v /i||£ 2 + 7 (VA V„/i)] 
satisfies (remember that a > j3) 

A\\h\\l 2 + ((3/2) [\\V x h\\h + W^Ml*] < m 

<A\\h\\l 2 + (3a/2) [||V^|| 2 + ||V^||| 2 ] . 

Moreover, since Hi(h) has zero mean on the torus, we have by Poincare's inequality: 

\\h\\l < 2 \\h - Tk(h)\\l + 2 WTUml < C - n,(fc)||i + l/2||V^||i 2 ) 

for some explicit constant C > 0, and similarly 

\\v x h\\l < c {\\v x h - n,(v x /i)||i + i/2||v^||i 2 ) 

for some explicit constant C" > 0. Hence we deduce that 

j/it)<-K (||/C + ||v x >||i) 

for some explicit K > 0, and that F{t) is equivalent to the square of the H 1 norm 
of h. We define the norm H 1 by 

|| • || w i = |A|| • ||| 2 + a|| V. • ||| 2 + /3||V„ • ||i 2 + 7 (V.-, V,-) L2 j . 
This concludes the proof. 

3. Some generalizations 

3.1. Higher-order Sobolev spaces. In this section we show how to extend the 
previous method to higher-order Sobolev spaces. Let us first introduce some no- 
tations about multi-indices. For a multi-index j in N N , we shall denote by 
the value of the i-th coordinate of j (i = 1, . . . , N), and by \j\ the I 1 norm of the 
multi-index, that is \j\ = Yl!i=i c i{j)- We also denote by 5i the multi-index such 
that Ci(8i ) = for i ^ i and Q (5j ) = 1. Finally for two multi-indices j and I in 
we set the shorthand &[ = d/dvj djdx\. 
In order to treat the higher-order derivatives we shall strengthen assumptions HI 
and H2 into 
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HI'. We assume HI. Moreover we assume that for any k > 1, for any multi- 
indexes j and I such that k — \j\ + \l\ and \j\ > 1, we have 



(djA(h),dfh) L2 >4\\djhf A -4\\h\\l 



2 



for some constants v$,Vq > 0. 

H2'. We assume that K has a regularizing effect in the sense: for any k > 1, 
for any multi- indexes j and I such that k = \j\ + |/| and \j\ > 1, for any 
5 > 0, there is some explicit C(5) such that 

(3.1) (djK(h), djh) L 2 < C{5) \\h\\ 2 Hk ^ + 5 \\d}h\\ 2 L2 . 

Again these strengthened assumptions are satisfied by the physical models we 
discussed in the introduction, as we check in Section Now we can formulate the 
coercivity estimate on T and the consequence on its semi-group. 

Theorem 3.1. Let L be an operator on L 2 satisfying assumptions Hl'-H2'-H3 
and k be an integer number. Then T = L — v ■ V x generates a strongly continuous 
evolution semi-group e tT on H k , which satisfies 

||e tT (Id-ry||^ <C T exp[-r T t] 

for some explicit constants Ct, Tt > depending only on the constants appearing 
in Hl'-H2'-H3. More precisely 

Vh€H k , {Th,h) Hk <-C' T [ J2 \\di(h-U g (h))\\: 

\\j\ + \l\<k 

for some (explicit) Hilbert norm 7i k equivalent to the H k norm, and some explicit 
constant C' T > 0. 

Proof of Theorem \3.1\ We write H k for the homogeneous Sobolev semi- norm, i.e., 

\\h\\% k = £ mi*. 

\j\+\l\=k 

Again we pick h e H k r\N(T) ± C\Dom(T) and we observe that h will stay in N(T) ± 
for all times. For k = 1 the result is given by Theorem ll.il We proceed by induction 
on k. 

First note that since the equation commutes with x-derivatives we have for the 
purely a> derivatives the analog of equation (|2.2|) . namely 

(3.2) ^h\\ 2 L2 < ^XWdfh-U^Wl 
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For derivatives including some f-component (i.e., \j\ > 1), by means of H2' with 
5 < (1^5 Vq )/(2vi ) we have the following estimate 

(3.3) 

±\\d{hf < -2 Yl (dih,d>- s 5 ;h) L . + (2C(5) + 2 4) \\h\\l^-4\\dih\\i. 

i\ci(j)>0 

For all I with \l\ = k and > 0, we consider the mixed term 

^ t (dt s h,d?h) = -\\d?h\\l 2 + 2 (L (dfh - U^h)),^). 
By means of HI we obtain 

(3.4) ±{df*_ Si h,d?h) < -\\dfh\\h + C L rj\\d¥h - H(S{70IIa + C L V - 1 \\dtsMl 

For / with |/| = k and i such that q(/) > 0, we define the following combination 
of derivatives of order and 1 in v : 

(3.5) Qi ti := a \\%h\\b + P ll^llk + 7 

By adjusting the constants a,j3, 7 > and using Poincare's inequality in the same 
way as in the proof of Theorem 11.11 in Section |21 it is straightforward to obtain 

(3.6) C (\\d°h\\i 2 + Hazily < Qu < c (||5>||| 2 + ii^vuy 

for some explicit constants c,C > 0, and that the time derivative fulfills the following 
inequality 

(3-7) ±Qu < -K {\\dMl + WdtsMD + C ° Ml*-* 

for some explicit constants K, Cq > 0. 

Now we combine all the derivatives in the following way 

(3.8) £ (D^n^ + l^)- 2 "- 1 ' yj a „ 

|i|+|il=fc,b1>2 |/|=fc,i| Cl (0>o 
where z/ = > is the constant such that (by assumption HI) 

\\h\\h > "0 \\h\\ L 2- 

By (|3.6jl . is equivalent to the square of the homogeneous Sobolev norm H k . 
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To estimate the mixed terms in the right hand side of (|3.8|) coming from (J3.3j) . 
we write 

-2s 



E 2(^y S (dih,di^;h) L2 < 

\l\ + \j\=k,\l\=s 

E (fP'iiwi 2 * E (IP'ik+Mi^ 



\l\+\j\=k,\l\=s \l\ + \j\=k,\l\~- 

\ E ®)~*m\+\ E (f P" nOiii. 

\l\ + \j\=K\l\=s |J| + | 3 -[=Ai,|{|= a 

and we derive by combining ()3.2|) . ()3.3j) and ()3.4|) that the time derivative of ()3.8j) 
satisfies 



V^)<^ + iw--eC oV2 ' z 



b'l>2 



^0 



-2(fc-l) 



Eiiwii+ 2 Eii« 

|il=i [i|=* 



By using ()3.6|) we end up with 



jT k (t)<C+\\H\^-K- \ E II W 



b'l+l'l=fc 



for some explicit constants C + , EE > 0. Together with the induction assumption for 
J-\, . . . , J-k-i this concludes the proof of the step k by considering some combination 

ol F k . □ 

3.2. Weak external potential. If the particles are subject to an external force 
field, which is given as the gradient of a scalar potential V, the evolution equation 
on the distribution generalizes in the following way 

(3.9) dtf + v V x f - V X V ■ V v f = Q(f, /). 

We still consider for the spatial domain Q = T N . We assume that V = V(x) is C 2 
and e~ v G L\. For simplicity we restrict to collision operators for which the problem 
without potential admits a Maxwellian equilibrium in this section and the next one. 
Moreover in this subsection we shall restrict further to collisional models admitting 
only mass conservation as a conservation law, i.e., we restrict to the relaxation 
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model. The stationary solution is determined by mas conservation and given by 

„ -y Poo . A 

foo = e - — V r— M 

e Wl 1 

II II J-'x 

where is the total mass of / and A4 is the normalized Maxwellian 

e -M 2 /2 

MM 



(2tt) n / 2 ' 

We also assume that the external field is weak in the sense: 

(3.10) H^llc 2 ^) < £ 

for some e depending on the collision operator. 

Then we consider fluctuations around equilibrium of the form 

/ foo v/oo h 

and we compute the following linearized equation on h: 

(3.11) d t h + v- V x h-V x V -V v h = L(h) 

where L is the linearized operator associated with Q as before. We define the 
operator T by 

T = L - v ■ V x + V X V ■ V„. 

Let us assume that L satisfies assumptions H1-H2-H3 and that the kernel of L in 
L\ is given by SpanjA^ 1 ^} (this assumption is satisfied for the classical relaxation 
model). In this case the kernel of T in L xv is trivially given by Spanje - ^ 2 Ad 1 / 2 } . 

We only sketch the proof of the following result: 

Theorem 3.2. Under the assumptions H1-H2-H3 on L, there is Eq > such that 
for any V G C 1 satisfying ()3.10|) with e < 6q, the operator T above satisfies the 
conclusion of Theorem li.il If moreover conditions Hl'-H2'-H3 hold for L, and 
V G C k+1 , then the operator T satisfies the conclusion of Theorem \3.1\ 

Proof. First let us recall that the L? norm of h is decreasing using the non-positivity 
of L and the skew symmetry of — v ■ V x + V ' X V ■ S/ v in this Hilbert space. 

We will only show how to establish the bound on first-order derivatives. The 
generalization to higher-order is straightforward. 

Let us consider h G iV(T)- 1 n Dom(T) fl v . The time evolution for the L 2 norm 
of h and the L? norm of its gradient in v are unchanged. For the gradient in x, 
equation (|2.2|) is replaced by 

^||V,/i||i 2 < -2 A \\V x h - Ui(y x h)\\ 2 A + 2e\\V v h\\h\\V x h\\l?- 
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Finally for the mixed term we have to replace (|2.4j) by 

^(va v v h) < - \\v x h\\l* + c L v II v a /i - n,( v x h) \\l 

+ c L v - 1 l|v„/i||i + £||v„/i||i 2 . 

Hence following the same arguments as in Section |21 we find the following differ- 
ential inequality on the quadratic form T: 

< -K{\\h\\\ + \\V x , v h\\l) + ae{ || V^|| L2 ||V,/i|| i2 + \\V v h\\h) 

for some explicit constant a, K > 0. Therefore it concludes the proof for e > small 
enough. □ 

Remarks: 

1. It may be possible that for a spatial domain Q = M. N a modified version of 
this strategy could be applied, assuming additionally that V satisfies a log-Sobolev 
inequality on Q. 

2. This subsection about weak external fields illustrates the fact that our method 
is robust, since it is based on a priori estimates, which remain true up to a pertur- 
bation. 

3.3. Self-consistent potential. Let us consider a collisional kinetic model for par- 
ticles which interact through collisions and also through a self-consistent potential. 
In this subsection we exclude the semi-classical relaxation collision operators. For 
the potential we consider the physically most common case of Poisson interaction. 
More precisely 

dtf + v- V x f -eV x V- V,../ = Q(f, f) 



(3.12) 



A x V = p-p , / Vdx = 



N 



where the coupling is via the density p{t,x) = f f(t,x,v) dv. The equation models 
particles interacting by binary collision or by scattering a background medium in 
thermal equilibrium, and at the same time interacting via electrostatic forces in the 
case e = —1 or by gravitational attraction in the case e = +1. Existence of global 
classical solutions in the large to the Vlasov- Poisson system in the torus has been 
proven in [2]. For the system including various collision operators, the existence of 
classical solutions has been established very recently in the articles discussed in the 
introduction. 

We consider the previous equation in the torus x £ for e = — 1 (electrostatic 
interaction). It admits a global Maxwellian equilibrium /oo, determined by the 
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conservation laws (for a more general discussion of the possible equilibria, we refer 
to jUIj for instance). 

Then we consider the linearization f = foo + \[Jooh around this equilibrium. 
Discarding the bilinear terms yields 

dth + v V x h - (v ■ V X V) fU 2 = L(h) 



A X V=([ hfj'dv) 



where L is the linearized collision operator associated with Q. 
Now let us denote by T p the operator on L 2 defined by 

(3.13) T p = L(h) - v ■ V x h + (v ■ V x V(h)) fU 2 . 

Then (defining il g as before) we have the following theorem 

Theorem 3.3. Let L satisfy the assumptions H1-H2-H3. Then the conclusion 
of Theorem li.il still holds true for the operator T p defined in ()3.13|) . If moreover 
L satisfies assumptions Hl'-H2'-H3 ; then the operator T p obeys the conclusion of 
Theorem \3.1\ 

Proof of Theorem YS.tA The proof is almost exactly the same as to the one of Theo- 
rem ll.ll Therefore we shall only indicate the differences in the estimates. Essentially 
the norm has to be modified in order to take into account the interaction energy. For 
the L 2 norm, one has by integration by parts and using that L is mass conserving 

^ t (\\hf L , + \\V x Vf L ,) <-2\\\h-^m\' 



Similarly one has on the gradient in x 

l<i,j<N 



d 2 V 


2 


dxidxj 


LI J 



< -2AIIVJI- 11^/1)111. 



dt 

For the time evolution of the gradient in v one has the additional term 



-2 J (V X V ■ V v fU 2 ) h 



and for the mixed term, one has the additional terms 

d 2 V 



-2 / \ x VhfU 2 + C j (j2 



ViVj 



hfU 2 - 



Since it is straightforward that all these additional terms as well as the gradient of 
V can be controlled by the L 2 norm of h, the end of the proof is straightforward as 
in Theorem ll.il □ 
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Remarks: 

1. We need the interaction to be repulsive in order to have the right sign for the 
interaction energy. For gravitational self-interaction potentials, our method would 
work as well under an additional assumption of smallness of the potential, i.e., 
assuming e > and e < eo with eo depending on the collision operator in (j3.12j) (the 
proof is similar to the case of a weak external potential). 

2. As noticed in ^S], the non-linear term arising from a self-interaction Poisson 
potential can be controlled, in the energy estimates for H k with k such that E(k) > 
N/2, by 

C\\h\\ H u \\h(v) l ' 2 \\ 2 Hk , 
using Sobolev embeddings and the straightforward elliptic estimate in the torus: 

V/>1, \\V x V\\ m <C,||/|| h «-i. 

Therefore one can extend the construction of smooth solutions near equilibrium 
in Theorem 11.21 to the case when a self-consistent Poisson potential is added, as 
long as the coercivity norm A of the linearized problem is stronger than the norm 
II M^) 1 ^ 2 ^ 2 - This is true, for instance, for the Boltzmann equation for hard spheres 
(considered in [18J), and also for the Landau equation with 7 > — 1 (see Subsec- 
tion E3J). 



4. Application to full non-linear models near equilibrium 

In this section we prove existence and uniqueness of smooth global solutions near 
equilibrium thanks to the coercivity estimates on the linearized models. This yields 
also explicit exponential rate of convergence to equilibrium. Obviously there is 
nothing else to prove for linear models (such as the classical relaxation or linear 
Fokker-Planck equation). Therefore let us assume that the collision operator is 
bilinear and let us denote the remaining term in the linearization process / = 
foo ~\~ fco h: 

T{h,h) = f^ 2 Q{f 00 h,f 00 h). 

In this section we shall consider a linearized collision operator L satisfying as- 
sumptions Hl'-H2'-H3 (we make the additional assumptions Hl'-H2' in order to 
get coercivity estimates in higher-order Sobolev spaces). Moreover we shall assume 
on the bilinear form V: 

H4. There is ko e N and Cp > such that for k > k , 

\\T(h,h)\\ Hk <c r \\h\\ Hk I J2 IIWI* 

\|j|+|2|<fc 
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Then we have 

Theorem 4.1. Let Q a (bilinear) collision operator such that 

(i) equation admits an equilibrium < G L 1 (T N x WL N ); 

(ii) the linearized collision operator L = L(h) around with the scaling f = 

foe + fU 2 h satisfies Hl'-H2'-H3; 

(iii) the bilinear remaining term V = T(h, h) in the linearization satisfies H4. 

Then for any k > k ( where k is defined in H4 ), there is e > such that for 
any distribution < fo G L 1 with 

\\Uo-foo)f^ l2 \\m<e Q 

there exists a unique global smooth solution < / = f(t,x,v) to equation (jl.lj) . 
which satisfies 

\\(ft-foo)f- 1/2 \\m<C e e-^ 

for some explicit constant Cq,Eq,t > 0, depending only on the constants appearing 
in Hl'-H2'-H3-H4. 

Proof of Theorem \4-l\ We explain the proof by a priori arguments, on a given 
smooth solution. The construction of positive solutions thanks to the estimates 
above is based on, by now standard, fixed point arguments (we refer the reader to 
[IH1IIH] for instance). 
The function h — (/ — f^) foo~^ 2 satisfies H g (h) = and it solves 

d t h = T(h) + T(h, h). 

Then we estimate the time evolution of the 7i k norm, defined in Theorem 13.11 

d 



Vj u .^\ 2 nk = 2(Th,h) nk + 2(T(h,h),h) Hk . 



We deduce that 



k <-C T [ £ \\djh\\l\+C £ \\T(h,h)\\l k + e\\h\\l k . 

j\+\l\<k 



Hence, by taking e small enough, 



< ] ,2 ^ Cr 

,b'l+l'l<* / Vb'l+I'l< fe 



This concludes the proof by maximum principle since the A norm controls the L 2 
norm. □ 
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5. Proof of the general assumptions for physical models 
5.1. Linear relaxation. We consider the linear relaxation equation in the torus 

1 



(5.1) 



dtf + v-VJ 



f(t,x,v m )dv* M(v) - f 



for x G and v G M. N (N > 1). Here k > denotes the Knudsen number and M. 
denotes the normalized Maxwellian: 



M(v) 



\v\ 2 /2 



(2tt) n / 2 

with mass 1, momentum and temperature 1. This equation preserves the total 
mass of the distribution 

Vt > 0, / /(£, x, v) dx dv = / f (x,v) dxdv 

Jj N xR N Jj N xR N 

but admits no other conservation law. For a given initial datum / > 0, it admits a 
unique global equilibrium = p^Ai, where is the total mass of / , defined by 



Pc 



T N xR N 



fo(x, v) dx dv. 



Finally let us add that the Cauchy theory is straightforward for (|5.1j) since it is 
linear (see |3] for instance for more details on this equation). 
We rescale the equation as 

/ = foo + ^/Uh = pM+ p 1/2 Mh, 

where M := \fM.. The equation for h reads 



d t h + v ■ V x h = — 



K 



ti M'dv')M — h 



-: L{h) 



where we have used the classical notation h' = h(v'). We split the operator L into 



L — K — A, K(h) = 

Therefore L satisfies HI taking | 
forwardly (with C(S) = 0) from 



-([ h'M'dv')M, A(h) = K- 1 h. 



\ L 2 , and assumption H2 follows straight- 



V v K{h) = - ( I h'M'dv') V V M. 
K \Jr n J 



Observe also that the strengthened assumptions HI' and H2', that are necessary 
to ensure decay in higher-order Sobolev norms, are satisfied straightforwardly. 
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The operator L is local in x and t. When x is fixed, it is well-defined and bounded 
on Ll, and it is self-adjoint non-positive on this space. More precisely its Dirichlet 
form is given by 



Therefore its kernel is N(L) = Span{M}, and we define 11; the (orthogonal) pro- 
jection on this space in L? v : 



Hence assumption H3 is satisfied. 

5.2. Semi-classical relaxation. Let us now modify the previous relaxation equa- 
tion in order to take into account quantum effects associated with the Pauli's ex- 
clusion principle. We thus consider the following semi-classical model (weakly non- 
linear) for charged particles (here e G {—1, 0, 1}): 



for x G T and v G R , where M. is the normalized Maxwellian as before, and we 
have used the shorthand /' = f(v'). Regardless of the choice of e this scattering 
operator is mass preserving but admits no other conservation law. The case e = 
is the standard linear relaxation model of the previous subsection. 

In the case e = 1 this operator is probably the simplest model describing a gas 
of fermions relaxing towards the thermodynamic equilibrium for a perfect fermigas, 
that is the Fermi-Dirac distribution. The operator Q+\ describes the interaction of 
the fermions with a background medium at rest with constant temperature. The 
factors (1 — /) correspond to correlation of particles before and after collision due to 
Pauli's exclusion principle. Note that this modification of the standard relaxation 
mechanism is (at least for the usual range of temperatures and densities) necessary 
only in particular situations, such as the one of a gas of electrons. Those are, 
due to their small mass, most likely to satisfy Sommerfeld's degeneracy condition 
(see |8j). This equation can be seen as the scattering counterpart of the full fermion 
Boltzmann equation studied for example in [HJEI]. For a more detailed introduction 
to models describing scattering as well as binary collisions for fermions see jHHl ■ For 
this model, a Cauchy theory can be obtained using maximum principle arguments to 





Then it is straightforward that L has a spectral gap Xl = k , since 




(5.2) d t f + vV x f 
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treat the (weak) nonlinearity (see [33]) assuming some bounds on the initial datum. 
The long-time behaviour of solutions to this equation has been studied by the EEP- 
method in [33] (leading to polynomial rates of convergence to equilibrium), however 
the necessary uniform regularity bounds on the solution were assumed. 

In the case e = —1, Q_i describes the interaction of bosons with background 
medium at rest with some constant temperature, and is probably the simplest model 
describing a gas of bosons relaxing towards the thermodynamic equilibrium for a 
perfect bosongas, that is the Bose-Einstein distribution. In this setting the existence 
of a boson in a velocity and space interval will increase the chance of another boson 
being scattered to this interval (see 8j). This mechanism leads to Bose-Einstein 
condensation for low temperatures and large densities. However since we linearize 
around the regular equilibrium we cannot describe this phenomenon but only the 
situation farther away from the critical mass for the phase transition. In the spa- 
tially homogeneous setting very precise asymptotics including optimal rates for the 
convergence above as well as below the critical mass have been given in ^3]. The 
authors study a model for Compton scattering of photons against electrons. In the 
non condensate case their model - in an appropriate scaling - is analogous to the 
one we study here because in this case their cross-section is bounded away from 
0. The authors prove exponential convergence in the non-condensate case, which is 
consistent with the result we obtain in the x-dependent situation. Moreover they 
derive an optimal rate for the convergence in the condensate case which is only poly- 
nomial. Thus is seems unavoidable to impose a bound on the initial mass to retain 
exponential convergence. The existence of solutions to a more elaborate collisional 
model (in the spatially homogeneous case) has been shown in [2E] (see also p9j). In 
the same work the weak convergence to regular equilibrium states (similar to the 
ones that our simplified model admits) has been shown to hold true above a specific 
temperature, larger than the critical one, and for which an explicit bound is given. 

The equation preserves the total mass of the distribution 



and admits (recalling in the boson case the bound imposed on the initial datum in 
Theorem II .3|) a unique equilibrium 



where is determined by the mass p of fo- 

We linearize the equation for a general scaling / = foo + m h (m is a given positive 
function only depending on v). Discarding the bilinear term, the equation for h reads 




f, 



oo 



1 + en^M ' 
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after straightforward computations 
1 



d t h + v ■ V x h 



m'(l + en^M') 1 + e/tooAl , 
^ — irrr-M ~ h m-M 



dv' =: L m (h). 



We make the choice m — (1 + tn^M.) 1 M^fk~^ l) where M := \[JA. It yields 
L = K - A with 



= - ( f h'M'dv') M, 



and A is the multiplicative operator by v with 

m = \ (1 + «~M> (/ r „ Y^^T */) = ^ (1 + eOM). 

Therefore L satisfies HI taking || ■ ||a = || • \\ L 2 , and assumption H2 follows straight- 
forwardly (with C(5) =0) from 

V v K(h) = - ( I h'M' dv') V V M. 

Again the strengthened assumptions HI' and H2' are also satisfied straightfor- 
wardly 

The resulting operator L is local in x and t. When x is fixed, it is well-defined 
and bounded on L\, and it is self-adjoint non-positive on this space. More precisely 
its Dirichlet form is given by 



(L(h),h)„ = -±- [ 



R N xR N 



h^ + en^M) h'(l + en^M') 



2 

X 



M M> 

(1 + en^M)- 1 (1 + tn^M'Y 1 M M' dv dv' 



-- f 



hM h'M'^ 1 



foo f'oc 



foe fL dv dv'. 



Therefore its kernel is 



N(L) = Span t % 



We define II; the (orthogonal) projection on this space in L? v : 

M ' 



n '(^(L^ d "') 



First let us explain how to show assumption H3 by non-constructive approach. 
As v > y_ with y_ = inf M jv v > 0, we deduce that the (bounded) multiplicative 
operator A on L\ has its spectrum included in (— 00, — y\. Then the operator K 
is straightforwardly compact on L\ and thus by Weyl's theorem about compact 
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perturbation of the essential spectrum of a self-adjoint operator in a Hilbert space, 
we deduce that the essential spectrum of L is included in (—00, —v\. The remaining 
discrete spectrum lies in R_ because of the sign of the Dirichlet form, and as is an 
isolated eigenvalue, we deduce that there is Ao > such that the non-zero part of 
the spectrum of L lies in (—00, — Ao]- This immediately shows H3 with Xl = \q/V 
where V = sup R jv v < +00. 

It is interesting to note that if the mass in the boson case matches the critical 
one then v is still nonnegative but it becomes zero for v = 0, which of course is the 
place where condensation happens. In this case our method breaks down because we 
lack the spectral gap (the essential spectrum reaches and one only expects some 
polynomial rates of convergence to equilibrium, which is consistent with |14j). 

Second we restrict to the fermionic case (e = 1) and we explain how to estimate 
explicitly Xl. Let us consider some function h orthogonal to M^ 1 f^. Then 



hM dv 



Hence we deduce that 



hMdv) < 



h M 



fo 



dv 



h 



I- 



,M 3 



1 + KocM 



dv . 



h vdv 



-1 "'oo 



dv 

h 2 v dv 



= 1 ( / p 

From the exact formula for it is straightforward that < and thus 

1 



fidv). 



P \Jr n 



fLdv) <1 



We get 



(L(h),h) 



L 2 



K \Jr n 



hM dv 



h vdv 



P 



< 



h 2 vdv I . 



Since in the fermionic case v > (kKqo) 1 , we deduce that L has a spectral gap Ao, 
with the explicit estimate: 



Xr, > 



fldv 



Hence we deduce an explicit estimate on A^ since A^ = Xq/v where v = sup R iv v 
(1 + K oo A^(0))/(/tfi; oo ) is explicit. 
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Now we want to establish the bound H4 on the bilinear part. It is given by 

K \J R N l + K^M' J 

Therefore H4 is immediately obtained by using Leibniz rule on higher-order deriva- 
tives, the trivial bound L 2 x L 2 — > L 2 on T, and Sobolev embeddings (which requires 
that E(ko/2) > N/2 where E denotes the entire part of a real number). 

Remark: The scaling that we used to linearize the collision operator is not exactly 
the same as in Theorem 14. II since we choose m = fli 2 (1 + e^oc^M) . However it 
is easy to see (following exactly the same proof) that the statement of Theorem 14.11 

— 1/2 

remains true also with the scaling / = f^ + m h when the factors foe are replaced 
by the some factors m _1 with the same decay at large velocities. This leads to the 
statement of Theorem 11.31 



5.3. The linear Fokker-Planck equation. We consider the linear Fokker-Planck 
equation in the torus 

(5.3) d t f + v V x f = V v ■ (V v f + fv) , 

for x G and v G R N (N > 1). 

This equation preserves the total mass of the distribution 

Vt > 0, / f(t,x,v)dxdv— / fo(x, v) dx dv 

Jj N xR N Jj N xR N 

but admits no other conservation law. For a given initial datum /o > 0, it admits a 
unique global equilibrium foo = p^Ai, where is the total mass of fo, defined by 

Poo = / fo(x,v) dxdv, 

JT n xR n 



and A4 is the normalized Maxwellian distribution 

MM ' 



\vf/2 



(2tt) n / 2 

with mass 1, mean and temperature 1. 

We study fluctuations around the equilibrium in the form 

/ = /«,+ y/JZh = Poo M + pV 2 Mh, 
where M := \^A4. The equation for h reads 



' N \v\ 

(5.4) dth + v ■ V x h = A v h + [ — - ^ ] h=: L{h). 
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We split the operator L into 

L = K — A, K = 0, A = L 

(K = is typical for a purely diffusive collisional model). Assumption H2 is obvi- 
ously fulfilled with C(5) = 0. 

Let us prove that L satisfies HI taking 

\\hh = {\\vh\\ 2 L2 + \\v v h\\h) 1/2 - 

Indeed one can check easily that this norm is stronger than the L 2 norm: 



\h\\i > 2 \\vh\ 



L 2 



\V v h\\ L 2 > -2 / (vh) -V v hdxdv = N\\h\\ 2 L2 . 



Straightforward computations yields 



(5.5) 



(A(h), h)ip. > Ci \\h\\\ — C 2 



|2 
\L 2 



for explicit constants C±, C 2 > 0. The operator L is local in x and t. When x is 
fixed, it is well-defined and bounded on L 2 ,, and it is self-adjoint non-positive on this 
space. More precisely its Dirichlet form is given by 



{L(h),h) Ll = - 



V v h + -h 



dv. 



Therefore its kernel is N(L) = Span{M}, and we define 11; the (orthogonal) pro- 
jection on this space in L 2 : 



ti M' dv' M. 



Classical computations based on Poincare's inequality with measure M. show that 



V v h + -h 



dv>2 \\h\\ 2 L 2. 



(5.6) 

Then combining ()5. 5I5.6|) yields 

(L(h), h) L 2 < —A \\h\\\ 

for some explicit constant A > 0. 
Finally we have 

(V v L(h),V v h) L 2 = (L(V v h),V v h) L 2 - ((v/2)h,V v h) L 2 

= (L(V v h), V v h) L 2 + — \h\ L 2 < -A L \\h\\\ 
The two last inequalities conclude the proof of HI and H3. 



N \\h\\ 

\\h \\t2. 

2 11 l|L 
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5.4. The Boltzmann equation. Let us consider the Boltzmann equation (here 
N>2) 

(5.7) dtf + v VJ = Q(f, /), t>0, xeT N , veR N 

with a collision operator (local in t, x) 

Q(fJ)= [ B(\v-v.\,cos9) (/'/:-//*) dv*da. 

We adopt the notations /' = f(v'), /* = /(i>*) and = /(i>*), where 

u' = (u + u*)/ 2 + (l u - v *l/ 2 ) a > < = (u + u*)/ 2 - (l u - v *\/ 2 ) a 
stand for the pre-collisional velocities of particles which after collision have velocities 
v and v*. Moreover 6 G [0, n] is the deviation angle between v' — v[ and v — v*, and B 
is the Boltzmann collision kernel determined by physics (related to the cross-section 
£(i> — i>*, cr) by the formula B — \v — v* \ E). On physical grounds, it is assumed that 
B > and I? is a function of |t> — i**] and cos0. 

Boltzmann's collision operator has the fundamental properties of conserving mass, 
momentum and energy 

/ Q(f,f)<j>(v)dv = 0, </>(v) = l,v,\v\ 
Jr n 

and satisfying celebrated Boltzmann's H theorem, which writes formally 

~f f\ gfdv = -[ Q(f,f)\og(f)dv>0. 
The equilibrium distribution is given by the Maxwellian distribution 

Poo ( "' 

M (Poo, Moo, Too) (v) = — — ^TTTTJ exp 



|2 



(277^)^ 2 

where p^, Woo, are the density, mean velocity and temperature of the gas 
Poo = / f(v)dxdv, Moo = — / vf(v)dxdv, 

Jt^xR^ Poo Jt^xR^ 

= -J— [ \uoo -v\ 2 f(v)dxdv, 

iv Poo Jt^xR^ 

which are determined by the mass, momentum and energy of the initial datum 
thanks to the conservation properties. 

The main physical case of application of this subsection is that of hard spheres in 
dimension N = 3, where (up to a normalization constant) 

(5.8) B(\v — i>*|, cos6>) = \v — v*\. 

More generally we shall make the following assumption on the collision kernel: 
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Bl. We assume that B takes the product form 

(5.9) B(|v-i;,|,cos0) = - 6(cos0), 

with $ and b non-negative and not identically equal to 0. This decoupling 
assumption is made for the sake of simplicity and could probably be relaxed 
at the price of technical complications. 

B2. Concerning the kinetic part, we assume $ to be given by 

(5.10) <$>(z) = C$ z'' 

with 7 e [0, 1]. It is customary in physics and mathematics to study the case 
when <&(v — v*) behaves like a power law \v — t>*| 7 , and one traditionally sep- 
arates between hard potentials (7 > 0), Maxwellian potentials (7 = 0), and 
soft potentials (7 < 0). We assume here that we deal with hard potentials 
(or Maxwell molecules). This assumption is crucial since for soft potentials 
with angular cutoff (see below), the linearized operator has no spectral gap. 



B3. Concerning the angular part, we assume that it is C 1 with the controls from 
above 

(5.11) Vz6[-l,l], b(z), b'(z)<C b . 

This implies in particular that B satisfies Grad's angular cutoff (see [16 ). 
Note that the smoothness assumption on b could be relaxed by using trun- 
cations and mollifications in the proof. 

When b is integrable on the sphere S N_1 (as here thanks to B3), we define 



| L 1( S JV-1) .- 



1*71 

I § N ~ 2 1 / 6(cos 9) sin N ~ 2 9d9<+oo. 
Jo 



Without loss of generality we set £b = 1 in the sequel. Then one can split the 
collision operator in the following way 

Q(g,f) = Q + (gJ)-Q-(gJ) 



Q + (gJ) = / $(\v-v*\)b(caa9)g'J'dv.&T. 

Q-(gJ) = / mv-v*\)b(coa9)g*fdv*da=(&*g)f. 
We introduce the so-called collision frequency 



(5.12) u(v) = / mv-v*\)b(cos6)M(v*)dv*d(T = (®*M)(v) 

and denote by vq > the minimum value of v. 
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Using the notation M = Ai 1 ^ 2 the linearized collision operator is given by 
L(h) = M' 1 [Q(Mh, M) + Q(M, Mh)} 



M I <$>{\v-v*\)b{cos6)M* 



K_ v_ 

ML + M' 



h* h 



df* da. 



L is self-adjoint on the space L v . It splits between a multiplicative part and a 
non-local part as follows 

L(h) = K(h) - A(h) with A(h) = u(v) h 

and 

K(h) = L + (h) - L*(h) with L*(h) = M [(hM) * $] 

and 

L + {h) = [ <$>(\v-v*\)b{cose) [h'Ml + KM'} M*dv*da. 

K is bounded and compact in L 2 ,, as proved in [T7j . 

From the classical spectral theory of L it is well-known that with the usual changes 
of variables 



(Lh, h) Ll = - 



4 



<&(|u - b(cos9) 



K_ h!_ 

Ml + W 



M M* dvdv* da < 0. 



This implies that the spectrum of L in is included in R_. Moreover the null 
space of L is 

(5.13) N(L) = Span {M, v x M, v N M, \v\ 2 M} . 
Using the fact that 

^ A (i + M) 7 <K^) <^ A (i + l^l) 7 

for some explicit constants z/ A , v% > 0, and that V v v E L°° with explicit bound 
since 7 G [0, 1], we deduce that assumption HI is satisfied with the norm 

\\h\\ A = \\h(i + \v\y/ 2 \\ L 2. 

Assumption HI' can be proved by similar arguments. 
Now we want to prove that 

(5.14) Vh±N(L), —(h, Lh) L 2 > X \\h\\\. 

Controls from below on the collision kernel are necessary to ensure the existence 
of a spectral gap for the linearized operator. Under our assumptions, the non- 
constructive proof of Grad shows that L has a spectral gap. Moreover explicit 
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estimates on the spectral gap have recently been obtained in jT] and extended to 
explicit estimates of the form (J5.14)) in |32]- Following these results L satisfies H3 
(for the norm A) with explicit bound. 

Now we fix some 5 > and check that L satisfies assumption H2 (assumption 
H2' can be proved by similar arguments). 

Concerning the part L*, this accounts essentially to Young's inequality. We easily 
compute the kernel of the operator 

L*h{v)= I h(v + V) k*(v, V) dV 
Jvm N 

with 

k*(v, V) = M{v) 1/2 $(|V|) M(v + V) 1/2 . 
We introduce the splitting k* = k*' s + k* ,r , with 

k* e > a (v,V)=l m > E} k*(v,V) 

where I denotes some mollified indicator function. This induces the corresponding 
decomposition L* = L* ,s + L*' r . It is straightforward that 

\\L*/\\ L ^ L 2 -^0 

and 

|V^' S | , \V v k* £ ' s \ , \V v k*/\ < C{e) M{Vf/\ 

Hence we deduce 

(5.15) \\V v L^ s h\\ L2 < C(e) \\h\\ L 2 
and 

(5.16) \\V v L*/h\\ L2 < 5 \\h\\ H i + C(e) \\h\\ L 2 
if e is small enough. 

Now we turn to the part L + . We follow Grad computations f7| Sections 2 and 3] 
(recalled also in Chapter 7, Section 2]) to compute the kernel of L + , and apply 
the same kind of estimates as in [3T]. We make the changes the variables 

• a E § N -\ v*eR N — ► u = (i/ - v)/\v' -v\e S^ -1 , E R N : the jacobian 
amounts to change b into 

1(9) = 2 N ~ l sm N - 2 9/2b(9); 

• then u E S^ -1 , E R N — > uo E S^" 1 , u = v - E R N : the jacobian is 
equal to 1; 

• then keeping u fixed, decompose orthogonally u = u u + W with u E K 
and W E ui^: the jacobian is equal to 1; 

• finally keeping W E V ± fixed, u E § N ~\ u E R — ► V = Uquj E R n : the 
jacobian is (1/2) |\^|-C^v-i). 
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Thus we get 

L + h(v)= [ h(v + V)k + (v,V)dV 

with 

~\W\ 2 - IVI 



k+{v,V) = cst\V\- {N - l) / $(^1^12 + \W\ 2 )b 



\W\ 2 + \V\ 2 , 
M(v + W) l/2 M(v + V + W) 1/2 dW. 

This kernel can be written as 



k + (v,V) = est M(Vy^\V\ 



l/4| T/ ,-(JV- 



- 1 ) f $(y/\V\ 2 + \W\ 2 )b( 
Jwev 1 - V 



\w\ 


2 _ 


y\ 


2 


\w\ 


2 + \ 


v\ 


2 



M(v + V + W/2) dW. 

Moreover it is shown in Chapter 7, Section 2] that 

Hi r + ll -R^oo n 

||J-H>ij-^ ||l 2 ^l 2 ► U. 

We use this to perform the splitting 

k+ = kp s + kp r 

with 

kp s (v, V) = X {H < £ - 1} X m > e} k + (v, V), 

where X denotes some mollified indicator function. The corresponding decomposi- 
tion of L is denoted by 



L+ = Lp s + Lp r . 



It is straightforward that 



and 



\V v kp s \, \V v kt' s \, \VvK' r \ <C{e)M{Vf/\ 

Hence we deduce 

(5.17) \\V v Lp s h\\ L2 < C{e) \\h\\v 
and 

(5.18) \\V v Lp r h\\ L2 < 5 \\h\\ m + C(e) \\h\\& 
as long as e is small enough. 

This concludes the proof by gathering (|5.15|) . ()5.16|) . ()5.17|) and ()5.18|) . 
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Finally let us consider the bilinear part given by 
T(hi, h 2 ) = M~ l [Q(Mh, Mh) + Q(Mh, Mh)] 

$(\v - v m \) b(cos9) M* [{hYM' - (/ii)*(/i 2 )] dv* da 

$(M)6(cos0)Af, [(hYM' ~ (hi)*(h 2 )} duda 



with the notation u = v — u*. We estimate 
V(hx, h 2 )(pdv 

<cy\\ L 2 ( / / KhYWihYJMiv + ^iuyduda 

R N xS N-l 



dv 



<C\\<p\\v / KhYfKhzYfivydvduda 

< C \\<p\\ L 2 (\\hWL2 \\h 2 \\ Av + \\hx\\ Av \\h 2 \\ L 2) , 

which implies 

l|r(/ii,/i 2 )IU2 < c (\\hx\\ L 2 \\h 2 \\ Av + ||/ii|| Al) WfaWiz) ■ 

Together with Leibnitz formula to differentiate T according to v and x and Sobolev 
embeddings this concludes the proof of H4 for E(ko/2) > N/2. 

5.5. The Landau equation. This subsection deals with the Landau equation (for 
N >2) 

(5.19) dj + v VJ = Q(f, /), t > 0, x e T N , v e R N 

which features the collision operator 

Q(f, f)(v) = V v ■ ( [ A(v- v„) [U (V/) - / (V/)J dv, 

where A(z) = \z\ 2 <&(\z\) P{z), $ is a non-negative function, and P(z) is the orthog- 
onal projection onto z L , i.e., 

(PC*))* = 5 ^ - j$ 

We use again the notation /* = /(i>*). This operator is used for instance in models 
for plasma. In this case the interaction among the particles is via the Coulomb 
potential and $(|z|) = \z\~ 3 in dimension 3. For more details see [J21 Chapter 1, 
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Section 1.7] and the references therein. Indeed in this case the Boltzmann collision 
operator does not make sense anymore (see jHj Annex I, Appendix]). 

Landau's collision operator has the fundamental properties of conserving mass, 
momentum and energy 



g(/,/)0(v)d« = O, <j>(v) = l,v,\v\* 
and satisfying Boltzmann's H theorem, which writes formally 

~f f\ogfdv = -[ Q(fj)\og(f)dv>0. 
The equilibrium distribution is given by the Maxwellian distribution 

■M(/P 00 ,'U 00 ,T 00 )(t;) = , ^ \ N lo ex P ^ ' 



where p^, Uoo, are determined as in the Boltzmann case. 
We make the following assumption on the collision kernel: 
LI. We assume <3> to be given by 

(5.20) $(z) =C< i ,z~ i 

with 7 G [—2, 1]. By analogy with the Boltzmann equation, one could say 
that this assumption covers hard and moderately soft potentials. 

We consider fluctuations around equilibrium of the form / = M + Mh. The 
linearized collision operator is given by 

'V v h s 



M 



MM^dv, 



L is self-adjoint on the space L\. It splits between an (almost) "convolution part" 
and a diffusive part: 

L(h) = K(h) - A(h) 

with 



K(h) = -M^Vv ( / A(v-«»)(^] MM-, dr.. 



and 



A(h) = -M~ l Vv^J n AO - v,) MM, dv, 

Estimate H2 on K is easily verified since 

K{h) = / h* dv* 
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where the kernel 



k(v,v*) 



(V„) 



belongs straightforwardly to L 2 



t ( M v - v*)MMi 
MM, 

and also to H 1 



(v„; 



X . 



X 



except possibly 



for a small region (v — v#) ~ which can split as in the Boltzmann case. 

It is well-known from the classical spectral theory of L that with the usual changes 
of variables we have 



(h,Lh) 



Li 



W - v* 



R^xR^ 



M 



M 



M M*dv dv* < 0. 
. Moreover the null 



This implies that the spectrum of L in l? v is included in K 
space of L is 

(5.21) N(L) = Span {M, Vl M, v N M } \v\ 2 M) . 

Now let us use some estimates proved in [T^j. First we define the norm 



h(v) 



1+7/2 



2 

hi 



(P(v)V v h)(v) 



1+7/2 



2 

hi 



{(l-P(v))V v h)(v) 



7/2 



hi 



2 



which is stronger than L 2 V as soon as 7 > —2. In T7B Section 2] it is proven that 

(Ah,h) L 2>C\ 

with explicit constant, from which HI follows. For the bilinear term Theorem 3 
from ^1 together with Sobolev embeddings yields H4 in the norm A with explicit 
constant as long as E(k /2) > N/2. Again the stronger assumptions Hl'-H2' 
are deduced straightforwardly with the same arguments. Finally assumption H3 is 
proved in |191 Section 2, Lemma 5] by non-constructive arguments (and an explicit 
proof is given in [33j). This concludes the proof. 

5.6. Remarks on other models. Linear models of radiative transfer in the torus 
enter straightforwardly our abstract framework. It is likely that linear scattering 
Boltzmann models or semi-conductors collisional models also do so. Moreover it 
is easy to see on the linear relaxation models (as well as on more general linear 
scattering models) that one could add with very few changes in our proof some 
scattering rate £ = E(x) depending on x in front of the collision operator: assuming 
that E G C°° and 

Vi E T N , < £_ < E(x) < £+ < +00 
for some constants £_, E + > 0, the conclusion of Theorem 11.11 still holds. 
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For the Boltzmann equation with soft potentials and Grad's angular cutoff, smooth 
solutions near the Maxwellian equilibrium have been built in t 20. : by including poly- 
nomial weight in v depending on the order of the derivatives in the energy estimates, 
it is likely that one can adapt our proof to build a norm which is decreasing along 
the flow inspiring from 20J. However in this case the integro-differential operator T 
is not coercive for this norm, instead it satisfies degenerated coercivity estimates for 
some weaker norms. This is enough to built smooth solution, but does not yield ex- 
ponential convergence towards equilibrium. Nevertheless as noticed in [37] one can 
deduce from it polynomial rates of decay to equilibrium by interpolating between a 
ladder of norms. 

Our analysis works at the linear level for the linearized Boltzmann equation for 
hard potentials without Grad's angular cutoff assumption, using explicit spectral 
gap estimates on L provided by [I] (note also that it could cover some moderately 
soft potentials interactions without Grad's angular cutoff, using explicit spectral 
gap estimates on L provided by jHH])- Indeed in [32| (see also jHHl), it is shown how 
to write the linearized collision operator in the form K — A with some regularizing 
K and some coercive A, and how to obtain coercivity estimates on L. However at 
now it is not known how to control the non-linear term in terms of a coercivity 
norm A adapted to the linearized operator. The functional space of these coercivity 
estimates in [32| is a local Sobolev space with the right fractional order, but which 
does not seem sufficient to control the non-linear term. 
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